HIP-1999-65/TH 

Quantum Magnetic Collapse 



M. Chaichian "'^ S.S. Masood'='^ C. Montonen^ 
A. Perez Martinez'^, H. Perez Rojas"'^'*^ 



"f^^"' High Energy Division, Department of Physics, University of Helsinki, 

. * Helsinki Institute of Physics, 

! P.O. Box 9, FIN-OOOI4 University of Helsinki, Finland 

\^ \ ^ Physics Department, Quaid-i-Azam University, Islamabad, Pakistan 

ICIMAF, Calle E No. 309, IO4OO La Habana, Cuba 
^ International Centre for Theoretical Physics 
^ ! P.O. Box 586, Strada Costiera 11, 34IOO Trieste, Italy 



(N 

We study the thermodynamics of degenerate electron and charged vector 
0^ ' boson gases in very intense magnetic fields. In degenerate conditions of the 

electron gas, the pressure transverse to the magnetic field B may vanish, 



X 



leading to a transverse collapse. For VF-bosons an instability arises because 
the magnetization diverges at the critical field Be = My^, je. If the magnetic 



Q \ field is self-consistently maintained, the maximum value it can take is of the 

' order of 2i?c/3, but in any case the system becomes unstable and collapses. 



Large magnetic fields can be generated due to gravitational and rotational effects in stel- 
lar objects like supernovas and neutron stars, i.e., magnetic fields of order 10^°G and larger 
have been suggested to exist in the cores of neutron stars [|I|. The standard electroweak the- 
ory establishes a limit on the magnetic field, the critical upper bound for stable vacua being 



Bf. = M'yy/e ~ 1.06 ■ 10^'' G, coming from the ground state energy eog = \J M'^ — eB, 
which is imaginary for B > Be. Fields of order Be may have been created at the electroweak 
phase transition (see 0, 0). The galactic and intergalactic magnetic fields can be consid- 
ered as relics of such huge magnetic fields in the early Universe @]- 0. In astrophysics, also 
the critical field Be' = m^/e ~ 4.41 ■ 10^^ G is relevant. 

Nielsen, Olesen and Ambj0rn [|T^ showed that the vacuum possesses the properties 
of a ferromagnet or an antiscreening superconductor for B ~ Be. It thus seems relevant 
to study the electroweak medium in a strong magnetic field of the order of the critical 
magnetic fields. The implications of these results for astroparticle physics and cosmology 
are expected to be interesting. As in preceding papers (Refs. [|l^, fl^), we only consider the 
first generation of leptons and quarks for the sake of simplicity. Here we shall calculate the 
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magnetization due to the charged leptons and intermediate vector bosons in the standard 
model. 

The thermodynamic potential Q = —TlnZ involves contributions from leptons and 
quarks, which are considered to be in chemical equilibrium among themselves through the 
boson fields, described by equations among their chemical potentials like = yUi/+/^e+, 
/^dz, + fJ'W+ = fJ'UL, fi'e+,w+ +^J'e-,w- = 0. From the thermodynamical potential we will choose 
the electron and W sectors exhibiting interesting effects in the astrophysical and cosmological 
scenarios respectively in the presence of extremely strong magnetic fields [B ~ B^' and 

In the astrophysical scenario the electron-positron gas thermodynamics is of interest. In 
the cosmological context, we will be concerned especially with the sector. 

The one-loop thermodynamical potential per unit volume of the electron-positron sector 

is Qe = rise + ^Oe, whcrC 



eB 



E 



(Eq+fle)|3■^ 



Here the sum extends over all Landau quantum numbers and the degeneracy factor is a. 
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q — Y pi + + 2ei?n and (3 = T ^. For VT' s, we have 
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where again we sum over all Landau quantum numbers and the degeneracy factor is b„ 
3-6on, with eog = 



w 



eB, and eg = Jpl + + 2eB{n + I). 



The Euler-Heisenberg vacuum terms are, for the electron-positron field. 
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and for the charged gauge bosons, 

e^B^ 
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ow 
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which diverges at i? > Be, leading to a vacuum instability. 

The mean density of particles minus antiparticles (average charge divided by e) is given 
by Ne^w = —d^e,w/df^e,w- We assume that there is always a background charge of opposite 
sign, to preserve electrical neutrality. We have 



eB ^ 



(5) 



where = [exp{Eq =f fie)P + 1] ^. 
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In the degenerate limit one gets A^'e = ^ X^o'' '^ny l^^l — — 2eBn, where the integer 
= I[{fil-m')/2eB]. 



For W-s, 



w 



eB 
4^ 
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h nj, = [exp(eog =F /ivy)/^ - 1] S = [exp(e, =F yUw/)/? - 1] ^• 

The magnetization is given by the contribution of electrons and charged vector bosons. 



It depends on the density of particles ylus antiparticles, and it is, 

M.w,e = -dQw,e/dB 
where (calling Moe,ow = -dQoefiw/dB), 



B 47r2 



, eBn, , 
dp3-—{n2+n^) 



+ M 
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and in the degenerate limit |T^, 

e / I 

A^e = — - a„ yUg-i/yUg — w? — 2eBn — (m^ + AeBn) In ■ 



-\/m^ + 2eBn 




and 



87r2 
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dp3 
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It is now especially interesting to discuss the equation of state of the system. The 
total energy-momentum tensor, whose spatial diagonal components are the pressures along 
the coordinate axes, may be obtained by starting from the quantum statistical average 

A" — 6nuC |]T3. If C is the total Lagrangian, after 



T^i. =< %v >s where T^^ 



dC 



doing the statistical average, its place in the energy-momentum tensor is taken by Q (since 
n = -p-Hn < e/o''^^4/d^x£(x4,x) ^j^g present SU(2)xU(l) model, the only non-zero 

averaged components of the field tensor are those of the U(l) external magnetic field tensor 



Ffj^p and then. 



T^, = {Tdn/dT + fiedn/dfie + fiwdn/dfiw)KS,^ + AF^pF,pdn/dF^ - (n) 

For Fpp = 0, ( p!TD reproduces the usual zero field case []T3|. For the electrically charged 
particles, we obtain thus different equations of state for directions parallel and perpendicular 
to the magnetic field. 
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P3 = p± = - BM. 



(12) 



This anisotropy in the pressures P3,J9_l leads to a magnetostriction effect in the quantum 
magnetized gas of charged particles. If ( |Tl| ) is taken as the Maxwell stress tensor (classical 
case), TVI < and > Ps, which produces a flattening effect in white dwarfs and neutron 
stars models |T^, In the present quantum case, for diamagnetic media also M. < Q 



leading again to a flattening effect. But for positive magnetization, the transverse pressure 
exerted by the charged particles is smaller than the longitudinal one by the amount BM.. 
The extreme case is found for magnetic flelds, eB ^ T^, when the electrons are conflned to 
the Landau ground state n = 0. (In what follows we will ignore the vacuum contribution to 
electron-positron pressure and magnetization, which is justifled at the scale of densities and 
flelds considered below). We have fig = —BAi^. where. 



Me = ^\^ieJ^Ji-m'^ -m^^n (13) 



and yUg — Y (27r2A^g/e-B)2 + m?, being the electron density. As nl > ni^, the expression 
(0) is always positive the system behaves as paramagnetic or ferromagnetic. But one of 
the most important effects we have in this limit is that the transverse pressure can vanish, 

P± = - BMe = 0. (14) 

(This is the lower bound for the pressure. For fermions, the pressure cannot be negative). 
The effect ([T^) is of pure quantum origin and it is easy to understand since all electrons 
are conflned to the Landau ground state, and the quantum average of their transverse mo- 
mentum vanish. If we consider a white dwarf star in which the predominating contribution 
to the pressure is from the electron gas, the vanishing of p± means that the gravitational 
pressure (of order GM'^/R^ where R is the geometric average radius of the star) cannot be 
compensated and an instability appears leading to a transverse collapse, i.e., the resulting 
object (a neutron star or a black hole) would be ellipsoidal, in this case stretched along the 
direction of the magnetic fleld, as a cigar. It is interesting to flnd the critical conditions for 
the occurrence of this conflnement to the state n = 0, and in consequence, for the collapse. 
We have. 



= = -r^ - 4.75 X 10-2°-|, (15) 



and the condition I{x) < 1 might be satisfled in some astrophysical conditions. E.g., for 
Ne ~ 10^°, B = 3.36 X 10^^ G, it is enough that B > Be' to satisfy it. For densities of the 
order of neutron stars, where a background of electrons and protons exist, if A^e = 10^^, the 
previous condition, if valid, would lead to i? > lO^^G. 



The W population in the Landau ground state is signiflcant ii d = J — eB < T. 



In the degenerate limit, e.g. for yM^ + eB /T ^ 1, one can neglect the contribution from 
excited Landau states and by taking only the n = term in (|lOl), one can approximate the 
flrst two terms, since the main contribution to the integrals comes from very small momenta. 
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The first term is the diamagnetic contribution which vanishes as T — >^ 0. The third is the 
vacuum contribution, which is asymptotically 

Mow--^-^-"^HM'w/eB-l), 

whose most important term is the second one which contributes para- or ferromagnetically 
for B > M^/2e, having a logarithmic divergence as 5 — >■ Be- As the logarithm is negative 
for Bc/2 < B < Be, that term has a negative contribution to the transverse pressure of 
vacuum for fields in that interval. The first term of Aiow contributes diamagnetically. But 
for B —>■ Be the dominant term in (|T6[) is the second, which is also para- or ferromagnetic, 
having a stronger divergence (inverse square root) than the vacuum term. To have a more 
explicit form for (|TB|), one must write /iiy in terms of the charge density. When confined to 
the Landau ground state the charge density of the system may be approximated as (po is 
some characteristic momentum) 

j^^^eBT p dp, rpo dp, eBT 2/.^ 

2vr2 ^0 Jp2 + rf2 _ Jo Jp2 + d^ + ^^ 47r 



from which /x^ = d'^/[l + and 



w 



Mw = 




Atx^Atx^N^ + e252T2 Ar^d V e^B^T^ 

+ Mow- (18) 

Taking > 10^^, T ~ 10~^ ergs and B < Be, one can neglect the unity in the square 
root and contributions from the first diamagnetic term in (p!6|) and from M.^, and one is left 
with 

A.. . H^. (19) 

The most important consequence is that the contribution of this magnetization to the 
transverse pressure of the W gas would be negative (see (|TU|)), and if M.wB contributes 
more than the pressure of other species, (the partial pressure p, = Qw even decreases as 
B — > Be) an instability occurs since the total pressure would be negative. Thus, for stability 
(also to prevent W decay), we must assume some background able to keep the total pressure 
P± > 0. 

Some sort of Bose-Einstein condensation actually takes place |T^ for bosons. For small 
momentum and magnetic fields strong enough B ~ B^, the term 1/d dominates and the main 
contribution to the W propagator comes from the low momentum gauge bosons |T^P |. 



In the absence of a magnetic field, the quantum degeneracy of the l^-boson sector leads to 
condensation, which at T ~ has been estimated [|T^ to occur induced by neutrino densities 
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of order Ucv = M^^/Gn"^ ~ lO^^cm"^. At any temperature, a spontaneous magnetization 
would appear in the condensate of charged bosons, say W~^, even at zero external field 
H = B — AnAi = 0. This spontaneous magnetization could self-consistently maintain the 
microscopic field B = inj^e^w- 

Let us assume the magnetization large enough to maintain the internal field B self- 
consistently and very large densities, such that fie ^ iti. The dominant term in ( p!oD is (|19D 
~ eNwfid. At B Be Gauss, we obtain that the self consistent critical field is reached at 
an electron density ~ 10^^ electrons/cm^. 

At such field intensities J^w diverges, but if we write the self-consistency condition for 
the W sector, we have 

B = AtiM = (20) 

(Jj 

Let us write eB = x^M"^ and since < x < 1, we easily get 

.^VT^^^'-^ = A. (21) 

As M^r/e^ ~ 10^^ cm~'^, even for Ny/ exceeding largely the nuclear density, A can be 
extremely small (For A ~ 1, Ny/ ~ 10^^ cm~'^. The horizon of events is ~ 5.6 cm.) By 
writing y = x^\J\ — x"^ we have a curve having an increasing branch starting from x = 
up to a maximum at xm = ^^2/3, y = Ai = 2^/3/9, compatible with a density Nw ~ 10^^ 
cm~ . We have also a decreasing branch from X — Xm 

to a: = 1, compatible with densities 
smaller than 10^^ cm~^. Thus, eq.(pT| will not have real solutions for A > Ai. However, for 



A < Ai we have two real positive solutions for x ( coinciding for = 2/3). For A <^ 1, 



these solutions are Xi = \jA + A? j2 and X2 = Vl — A"^. The first solution means that B 
increases with increasing Nw, (up to the value Bm = 2Miy/3e). In the second solution B 
decreases as a function of Nw, its limit for Nw being Be- This obviously indicates 
that the expression for the magnetization must include the contribution from Landau states 
other than the ground state, which leads to a diamagnetic response to the field. This would 
compensate the increase of the self-consistent field with increasing Nw to keep B < Be- 

This can be shown to occur from formula ([TD|). If we call the ground state density 
Nwg and the density in other Landau states Nwn {Nw = Nwg + Z^-^VKn), for B > Bm = 
2Bc/3, dB/dNwg < and dB/dNwn > and excited Landau states start to be populated. 
The condensate in the ground state decreases in favor of the increase of the population 
in excited Landau states, which starts to grow and contribute diamagnetically to the total 
magnetization keeping AttA4 = B < B^. But for the system to react in this way, an enormous 
amount of energy (and angular momentum) would be required, of the order of respectively 
NwMw and Nw (here we neglect the running of Mw)- But the transverse collapse takes 
place at such densities: since the pressure comes essentially from the fermion (electron) 
background, the self-consistency condition B = 4:7iJ\4e,w leads to ps = —fl ~ BMlf. and 
p± = Ps — B{J\4e + M-w) ~ and thus the system collapses. 

Let us assume that in some stage of the early universe a very large external field H Be 
was present. If T ~ Mw-, as happened near and below the electroweak phase transition, 
using up the energy and angular momentum of the background radiation, pairs will be 
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produced in the energetic more favourable Landau ground state (having a "mass" = d), and 
this process would be even more favoured as the magnetic field approaches Be even for lower 
temperatures. The magnetization Ai is given by an expression similar to the second term in 
(p!0|), in which the expressions for the particle-antiparticle densities would be in equilibrium 
with the electromagnetic background radiation. This means that one must take the chemical 
potential as zero (equal number of W^). Then, A4w ^ e^BT/in'^d. The density of particle 
+ antiparticle pairs would then be ~ eBT ~ 10^^ cm~^, and the microscopic field B < Be 
starts to be maintained self-consistently. We would have the situation discussed in the 
previous paragraph. The process of W pair creation in the external field would lead to a 
creation of order from disorder, i.e., to an effective cooling of the subsystem considered, 
although due to similar reasons as before, p± and the system becomes unstable and 
collapses. 

We conclude, first, that if a degenerate electron gas is confined to its Landau ground 
state, its transverse pressure vanishes. This phenomenon establishes a limit to the magnetic 
fields expected to be observable in white dwarfs, and even in neutron stars. Second, that the 
instability of the vacuum in magnetic fields i? ~ i?c in a hot and dense medium, is avoided, 
since the self-consistent magnetization prevents fields greater than 2Bc/3, although under 
such conditions the system becomes also unstable and collapses. 
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